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Abstract 



O 
U 

Oh! 

Q ■ We present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama 

^ , Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (AC A). In its most 

c/^ ■ compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular 

C^l resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons produced 

during cluster mergers at high redshifts, that are hard to spatially resolve by current and near- future X-ray 

^S| \ detectors. Quality of image reconstruction is poor with the 12m array alone but improved significantly 

^ . by adding AC A; expected sensitivity of the 12m array based on the thermal noise is not valid for the 

On ■ Sunyaev-Zel'dovich effect mapping unless accompanied by an ACA observation of at least equal duration. 

vN , The observations above 100 GHz will become excessively time-consuming owing to the narrower beam size 

and the higher system temperature. On the other hand, significant improvement of the observing efficiency 

is expected once Band 1 is implemented in the future. 

CO . Key words: cosmology: observations - galaxies: clusters: - radio continuum: galaxies - techniques: 

^^ ' interferometers 

• • 

_ ^_^ 1. Introduction growth of cosmic structures was faster than today. 

^ I As of today, however, spatial resolutions of the major- 

^ . The Sunyaev-Zel'dovich effect (SZE, Sunyaev & ity of the SZE images are at arcminute scales with a lim- 

_C^j Zel'dovich 1972), inverse Compton scattering of the cos- ited number of exceptions (Komatsu et al. 1999, 2001; 

mic microwave background off hot electrons, provides a Pointecouteau et al. 2001; Mason et al. 2010; Zemcov et 

unique probe of the intracluster medium (for reviews see al. 2010; Massardi et al. 2010; Korngut et al. 2011). This is 

Rephaeh 1995; Birkinshaw 1999; Carlstrom, Holder & largely due to the low surface brightness of the SZE and a 

Reese 2002). For the gas with given electron density n^ small number of radio telescopes with sub-arcminute res- 

and temperature Tq at redshift z^ surface brightness of the olution in the range of wavelengths relevant to the SZE 

thermal SZE is proportional to the line-of-sight integral measurements. Controlling systematics of such large tele- 

of electron density times temperature, JueTedl, whereas scopes to the sensitivity level required for the SZE imaging 

that of the X-ray thermal bremsstrahlung emission is pro- is also challenging. Interferometers are therefore a promis- 

portional to (1 + z)'^ JnlTe^'^dL The relative significance ing and complementary tool for high sensitivity SZE imag- 

of the SZE over the X-ray emission thus increases with ing observations (e.g., Jones et al. 1993; Carlstrom, Joy 

redshift and electron temperature. This makes the SZE & Grego 1996; AMI Cohaboration 2006; Muchovej et al. 

a powerful tool in detecting the shocks and the hot gas 2007; Wu et al. 2009). Since it measures synchronized 

associated with violent cluster mergers (Kitayama et al. and correlated signals among different telescopes, an in- 

2004; Ota et al. 2008; Korngut et al. 2011), the frequency terferometer has a much better control of systematic noise, 

of which is expected to increase at high redshifts when the e.g., from the atmosphere. The obtained data also span 
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a wide range of spatial scales, making it possible to sepa- 
rate point-like sources from the SZE. On the other hand, 
it often has a limited field-of-view of 0.5' ~ 2' in the mil- 
limeter/submillimeter bands and requires a large number 
of pointings (or mosaics) to uncover the bulk of the clus- 
ter emission that extends over several arcminutes. The 
sensitivity also varies with spatial scales depending on ob- 
serving conditions such as array configuration. Addition 
of single-dish data further improves the sensitivity for ex- 
tended sources (e.g., Vogel et al. 1984; Stanimirovic et al. 
1999; Takakuwa et al. 2003; Kurono, Morita & Kamazaki 
2009). It is hence not always easy to know a priori to what 
extent a specific interferometer can be used to image the 
SZE within a realistic observing time. 

In this paper, we investigate the feasibility of 
SZE imaging observations by the Atacama Large 
Millimeter /submillimeter Array (ALMA).^ ALMA wih 
consist of at least 66 antennas, designed to operate at 
photon frequencies between 30 GHz and 950 GHz. A ma- 
jor array of 12-meter antennas will be combined with the 
Atacama Compact Array (ACA, Iguchi et al. 2009) made 
up of 7- meter antennas and 12-meter single-dish antennas 
to improve the spatial frequency coverage and the total 
power measurement. While the inclusion of ACA should 
greatly enhance the capability of ALMA in observing ex- 
tended sources, the practical feasibility and how to op- 
timize it with such heterogeneous arrays in real observa- 
tions are not straightforward and certainly depend on the 
nature of the target. We thus perform detailed imaging 
simulations of the SZE observations with ALMA includ- 
ing ACA using their latest configurations (see Pety, Gueth 
& Guilloteau 2001; Heifer et al. 2002; Pfrommer, Ensslin 
& Sarazin 2005; Takakuwa et al. 2008 for earlier imaging 
simulations of different targets). 

We pay particular attention to resolving the shock 
structures of distant merging clusters that are relatively 
compact and hence suit the field-of-view of ALMA. We 
take the so-called bullet cluster, IE 0657-558 at z = 0.296, 
as a representative example and perform mock obser- 
vations using snap shots of state-of-the-art numerical 
simulations for this cluster (Takizawa 2005; Akahori & 
Yoshikawa 2012). IE 0657-558 is well known for its promi- 
nent bow shock observed by Chandra (Markevitch et al. 
2002) with the inferred Mach number of 3.0 ±0.4 and the 
shock velocity of ~ 4700 km s~^ (Markevitch & Vikhlinin 
2007). It has also been observed frequently via the SZE 
(e.g., Andreani et al. 1999; Halverson et al. 2009; Plagge et 
al. 2010; Zemcov et al. 2010; Malu et al. 2010). The tem- 
perature or the energy distribution of the shocked elec- 
trons, however, is still poorly constrained even for this 
well-observed cluster because of the lack of sensitivity of 
current X-ray and SZE detectors. We will show that 
ALMA is a powerful instrument capable of resolving a 
merger shock structure of galaxy clusters at high redshifts 
and improving our understanding of cluster evolution. 

Throughout the paper, we assume a standard set of 
cosmological parameters: ftm = 0.27, I^a = 0.73, and h = 
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Fig. 1. The compact array configuration of 12m x 50 and 
7mxl2 used in our mock observations. The projected po- 
sitions on the UTM coordinate are shown with the symbol 
sizes representing the diameters of the telescopes. 

0.70 (Komatsu et al. 2011). In this cosmology, an angular 
size of V^ corresponds to a physical size of 4.4 kpc and 8.1 
kpc at z = 0.3 and ^ = 1, respectively. 

2. Imaging Simulations 

The imaging simulations are performed using the 
MIRIAD software (Sault, Teuben & Wright 1995), ver- 
sion 4.2.3. Following the currently planned specifications 
of ALMA, we consider an array made up of fifty 12-meter 
antennas (hereafter 12mx50) and ACA made up of twelve 
7- meter antennas (hereafter 7mxl2) and four 12-meter 
single-dish antennas (hereafter SDx4). A summary of 
telescope specifications are given in Table 1. 

To maximize the sensitivity for the extended signals, 
the most compact array configurations, shown in Figure 
1, are adopted for 12mx50 and 7mxl2.^ The simulations 
are performed at 90GHz, which is the lowest frequency 
(and hence has the largest field-of-views and the lowest 
system temperature) in the initially planned specifications 
of ALMA. This frequency is also expected to be an opti- 
mal choice for the SZE observation owing to the minimal 
contamination by synchrotron and dust emissions. The 
prospects for the observations in the other bands will be 
discussed in §4. 

2.1. Input Models 

We use three input models whose observational set-ups 
are summarized in Table 2. For each of the model im- 
ages, we adopt the pixel size of 1^', to satisfy the Nyquist 
condition for the smallest spatial scale, ^ 3^' full width at 
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The antenna configuration files are identical to alma. out Ol.cfg 
and aca.i.cfg used in CAS A 3.3. 



No. 



Simulations of the SZE for ALMA 



Table 1. Summary of telescope parameters. 



Number of antennas 12mx50, 7mxl2, SDx4 

Configuration Most compact 

Frequency 90 GHz (Band 3) 

Bandwidth 7.5 GHz per polarization 

Intermediate frequency 6 GHz 

System temperature 73.5 K 

Aperture efficiency 0.71 

Correlator efficiency 0.88 

Primary beam FWHM 69'' (12mx50), 118'' (7mxl2) 

Phase error (rms) 20° 

Gain error (rms) 0.1% 

Pointing error (rms) 0.6" 



half maximum (FWHM), probed by the longest baseline 
of ALMA at 90GHz in the most compact configuration. 
Wherever necessary, linear interpolation is performed to 
assign the intensity to each pixel. The model images are 
created to cover the entire observing area. 

2.1.1. Model A: Gaussian 

The first model is a simple two dimensional Gaussian, 
placed at zenith (RA,DEC) = (00 : 00 : 00, -23 : 00 : 00), 
with FWHM ranging from 5" to 60". The peak intensity is 
fixed at — 100/iJy/arcsec^, corresponding to the Compton 
y-parameter of 5 x 10~^ at 90GHz. We use this model sim- 
ply to examine how the feasibility depends on the source 
size when the noises are negligible. 

2.1.2. Model B: Bullet cluster 

The second model is a latest result of 3D hydrodynamic 
simulations of a cluster merger by Akahori & Yoshikawa 
(2012) that are designed to reproduce the observed prop- 
erties of the bullet cluster, IE 0657-558. They carried 
out a set of N-body and SPH simulations of a collision of 
two galaxy clusters with virial masses of 1.5 x 10^^ and 
2.5 X 10^"^ M0 and the initial relative velocity of 3000 
km s~^, using a code developed by Akahori & Yoshikawa 
(2010). We use the result after 1.12 Gyr from the ini- 
tial condition for their two temperature model and com- 
pute the intensity map of the SZE, including relativis- 
tic corrections (Itoh & Nozawa 2004; Nozawa, Itoh & 
Kohyama 2005). To improve the agreement with the ex- 
isting data, we renormalize the electron density so that 
the peak y-parameter (oc tIqTq) matches the central value 
of 3.31 X 10-^ inferred by APEX-SZ (Halverson et al. 
2009). This is realized by reducing the electron density 
in each mesh by 26%.^ To cover both the SZ emission 
peak and the shock front, the image center is shifted by 
— 15" {QQ kpc) and +20" (88 kpc) along the right ascen- 
sion and declination, respectively, from that plotted in 
Figs 1-6 of Akahori & Yoshikawa (2012) and placed at 
(RA,DEC) = (6 : 58 : 27.4, -55 : 56 : 50). As the spatial 
resolution of the simulation is moderate near the shock 



^ Reducing only the electron density yields the most conservative 
estimate of the predicted SZE intensity for a given y-parameter, 
because the relativistic correction leads to the larger reduction 
in the predicted intensity for the higher electron temperature at 
90GHz. 



front (26 kpc or 6.0"), we use this model mainly to ex- 
amine the feasibility of observing the global structure of 
a real cluster. 

To further examine the feasibility of separating radio 
sources from the SZE, we add to Model B point sources in 
the field of IE 0657-558 reported in the literature (Liang 
et al. 2000; Wilson et al. 2008; Malu et al. 2010). The 
source flux at 90GHz is estimated by the following proce- 
dure; i) wherever available, we extrapolate the observed 
fluxes at the two nearest frequencies to 90GHz assuming a 
single power-law, ii) if the observed flux is available only at 
one frequency, we extrapolate it assuming a single power- 
law with a spectral index of —0.5, iii) to be conservative 
against uncertainties in the extrapolation, we multiply the 
extrapolated fluxes by a factor of 2. The resulting source 
properties are listed in Table 3. 
2.1.3. Model C: Shock front 

The third model is based on a higher spatial resolu- 
tion Eulerian mesh simulations of a moving substructure 
within a main cluster by Takizawa (2005). The model 
used in this paper is nearly the same as "the radial infall 
model" in Takizawa (2005) except that the gas density is 
reduced by a factor of 10 so that the y-parameter around 
the shock front matches that of Model B. We compute 
the SZE map from a snapshot at the elapsed time of 0.78 
Gyr from the start of the simulation, which roughly re- 
produces the observed morphology of IE 0657-558, and 
place it at (RA,DEC) = (6 : 58 : 17.3,-55 : 56 : 30). This 
model predicts the y-parameter gap across the shock front 
of 1.5 X 10~^, about half the central y-parameter of this 
cluster. As this simulation has a higher spatial resolu- 
tion (2 kpc or 0.5") and a smaller box size (0.8 Mpc or 
180") than Akahori & Yoshikawa (2012), we use Model C 
to examine the feasibility of resolving a merger shock in 
a cluster. As our mock observations are performed on a 
slightly larger area than the box size of Takizawa (2005), 
we pad mirror images to the original simulation output to 
cover the entire observing area. 

2.2. Mock visibilities 

Mock interferometric observations are performed over 
the sky area of 120" x 120" (Models A and C) or 240" x 
240" (Model B). Linear mosaicing (e.g., Sault, Staveley- 
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Parameter 




Model A 


Model B 


Model C 


RA (J2000) 




00 : 00 : 00 


6 : 58 : 27.4 


6:58:17.3 


DEC (J2000) 




-23 : 00 : 00 


-55 : 56 : 50 


-55 : 56 : 30 


Hour Angle 




-1 hr ~ 1 hr 


— 5 hr ~ 5 hr 


-5 hr ~ 5 hr 


Mapping area 










12mx50, 7mxl2 




120" X 120" 


240" X 240" 


120" X 120" 


SDx4 




190" X 190" 


310" X 310" 


190" X 190" 


Pixel size 










12mx50, 7mxl2 




1" 


1" 


1" 


SDx4 




10" 


10" 


10" 


Number of mosaics/pointing 


s 








12mx50 




19 


67 


19 


7mxl2 




7 


23 


7 


SDx4 




400 


1024 


400 


Total integration time 










12mx50 




2 hr 


10 hr 


10 hr 


7mxl2, SDx4 




2 hr 


40 hr 


40 hr 


Noise added 










12mx50, 7mxl2 




thermal 


thermal, phase 
gain, pointing 


thermal, phase 
gain, pointing 


SDx4 




thermal 


thermal 


thermal 


Synthesized beam (major, minor, p. a.)* 








12mx50 




4.0", 3.8", - 7.3° 


4.9", 4.2", - 88° 


4.8", 4.2", - 88° 


7mxl2 




17", 16", 5.1° 


19", 18", 83° 


19", 18", 84° 


12mx50+7mxl2 




4.1", 3.8", -7.3° 


5.2", 4.6", - 88° 


5.2", 4.5", - 88° 


SDx4 




69", 69", 0° 


69", 69", 0° 


69", 69", 0° 


Thermal noise (/xJy/arcsec^, 


/iJy/beam)^ 








12mx50 




1.1, 19 


0.62, 15 


0.35, 8.0 


7mxl2 




0.48, 150 


0.15, 60 


0.099, 38 


12mx50+7mxl2 




1.0, 18 


0.52, 14 


0.30, 7.8 


SDx4 




0.64, 3400 


0.23, 1200 


0.14, 770 



Three numbers for each model denote the major axis FWHM, the minor axis FWHM, and the position angle, respectively. 
^ Two numbers for each model denote the rms values in /iJy/arcsec^ and yuJy/beam, respectively. 



Smith & Brouw 1996) is performed with pointing centers 
lying on triangular grids separated by half the primary 
beam FWHM given in Table 1. Each mosaicing center is 
observed repeatedly for 30 sec on-source and 5 sec toward 
an off-source calibrator with a sampling interval of 5 sec. 
The numbers of mosaics, the hour angles, and the total 
integration times are listed in Table 2. The target eleva- 
tion is above 20° throughout the hour angles considered in 
this paper. As shown in Appendix 1, effective integration 
time toward each sky point apart from the map edge is 
nearly uniform and given by 

teff^2.6-^, (1) 

where tint is the total integration time, TVmos is the number 
of mosaics, and the numerical factor accounts for the over- 
lap of mosaics for the grid orientation mentioned above. 

The visibility data are then created separately for 
12mx50 and 7mxl2 by performing the 2D Fourier trans- 
formation to the u — v plane using the MIRIAD task uvgen 
and uvmodel The visibility loss due to shadowing by ad- 
jacent telescopes is also taken into account. The thermal 
noise is computed in all models using the same param- 



eters as the ALMA sensitivity calculator ^ as of March 
2012; the system temperature, the aperture efficiency, the 
correlator efficiency, the bandwidth, and the intermediate 
frequency are listed in Table 1. The system temperature 
is computed by adopting the source declination of Model 
B (— 55 : 56 : 50) and precipitable water vapour of 2.748 
mm at 90GHz. 

For Models B and C, antenna phase and gain noises 
with rms values of 20 degrees and 0.1%, respectively, are 
also added over an interval of 5 minutes using the MIRIAD 
task gperror. Relative pointing errors are also added to 
these models in an approximate manner by convolving the 
input images by a Gaussian with 0.6^' rms before creating 
the visibilities. This would be sufficient for our present 
purpose because the pointing errors practically affect our 
results only via identifications of point sources and the 
u — v coverage of ALMA is very good as shown below. 

Figure 2 illustrates the u — v coverage of our mock ob- 
servations. The spatial frequencies that correspond to the 
baseline lengths from 4 to 48 kA and 2 to 10 k A are covered 
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No. 



Simulations of the SZE for ALMA 

Table 3. Positions and fluxes of model point sources. 





Input 






Reconstructed 




AR.A.[arcsec 


ADec arcsec 


Flux [mJy 


AR. A. [arcsec 


ADec 


arcsec 


Flux [mJy 


83 


-12 


1.0 


83 




-12 


0.93 


89 


-35 


0.74 


89 




-35 


0.69 


-108 


147 


0.47 


-108 




147 


0.42* 


144 


86 


0.45 


144 




86 


0.43* 


117 


-111 


0.43 


117 




-111 


0.37 


-68 


-111 


0.37 


-68 




-111 


0.31 


78 


113 


0.16 


79 




113 


0.14 


-95 


-94 


0.13 


-95 




-94 


0.14 


-25 


97 


0.079 


-25 




97 


0.092 


-30 


9 


0.035 








undetected 


54 


-68 


0.032 








undetected 


1 


44 


0.028 








undetected 



Sources lying outside the 240^^ x 240^^ region mapped by 12m x 50. 



nearly uniformly by 12mx50 and 7mxl2, respectively, 
where A = 3.3 mm at 90GHz. Properties of the result- 
ing synthesized beams are given in Table 2. The effective 
area of the synthesized beam is given by 



^eff,SD '■ 



Hnt,SD 



Ae 



81n2 



: 28 arcsec 



Ohe 



b" 



(2) 



where ^beam is the geometric mean of the major axis 
FWHM and the minor axis FWHM. 

In addition, single-dish images are created in real space 
assuming a Gaussian beam with 69^' FWHM and tak- 
ing account of the thermal noise. An additional factor 
of \/2 is multiplied to the thermal noise level so as to 
take account of beam switching. Since the spatial reso- 
lution of the single-dish data is considerably worse than 
that of 12mx50 and 7mxl2, the single-dish mapping is 
performed over a larger sky area to reduce the edge ef- 
fect; an extra 35^' (0.5 FWHM) guard band is observed 
on each side of the image to cover 190^' x 190^' (Models A 
and C) or 310'' x 310'' (Model B). As the primary beams 
of 12mx50 and 7mxl2 also extend into the guard band, 
we will be able to detect point sources in this band to 
some extent as will be shown in §3.2. The pixel size of the 
single-dish image is taken to be 10" (Table 2). While this 
somewhat oversamples the data, we adopt this pixel size in 
the present analysis for the following reasons; i) as the ex- 
pected random noise of the single-dish image reduces lin- 
early with the pixel size, the above choice serves as a con- 
servative limit, ii) the thermal noise of the single-dish data 
in /iJy/arcsec^ is already smaller than that of 12mx50 
with the above pixel size (Table 2), and iii) we find that 
adopting larger pixel sizes often results in spurious distor- 
tions in the reconstructed images after deconvolved jointly 
with the interferometric data. Yen et al. (2011) also show 
that one may well reproduce the extended emission by 
adopting pixel sizes that are sufficiently smaller than the 
primary beam of an interferometer. Effective integration 
time at each point of the single-dish image, including the 
time spent for beam switching, is given by 



N. 



point, SD 



(3) 



where tint,SD is the total integration time and Appoint, sd is 
the number of pointings of SDx4 listed in Table 2. 

The mock data so obtained are reduced without a priori 
knowledge of the input model or point sources as described 
in the following subsections. The point source removal 
and image deconvolution are both performed over the ex- 
tended regions with the guard band mentioned above. 
Impacts of other sources of noise than considered here 
will be discussed in §3.4. 

2.3. Point source removal 

Whenever point sources are present, we make use of 
long baseline data to remove them before reconstructing 
the extended emission. First, a dirty map and a dirty 
beam are created only for 12mx50 by an inverse Fourier 
transformation of mock visibilities using the MIRIAD task 
invert 

Second, the position and the apparent flux of a source 
whose peak brightness lies above a certain threshold 
value are identifled with the CLEAN algorithm by Steer, 
Dewdney & Ito (1984) using the MIRIAD task mossdi. 
If there are multiple sources, the flux of each source may 
be mutually influenced via the side- lobe effects. A single 
point source may also be detected as a superposition of 
separate sources. In either case, the apparent ffux of the 
i-th source can be expressed as 

N 



Es^^Fr 



where N is the total number of detections. 



pint 



(4) 



is the 



intrinsic ffux of the j-th source, and the matrix Sij de- 
scribes the contribution of Fj"* to F^^^ for a given beam 
shape and source positions. Given knowledge of the u — v 
coverage of the observations, one can estimate the com- 
ponents of Sij and solve the linear algebraic equations for 
Fj"*. In the present analysis, all the detections that lie 
within the synthesized beam FWHM from each other are 
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regarded as components of a single source, whose flux and 
position are given by the sum and the flux- weighted mean, 
respectively, of fluxes and positions of all the components. 
Finally, contribution of identified sources is removed 
from the visibility data of both 12mx50 and 7mxl2 as 
well as from the SDx4 image. 

2.4' Image deconvolution 

Having removed the identified point sources, we cre- 
ate dirty images and dirty beams for either each of or 
the combinations of 12mx50 and 7mxl2. The Briggs 
robust weighting scheme (Briggs 1995) is adopted with 
the robustness parameter of 0.5. The dirty images are 
then deconvolved (corrected for the effects of the dirty 
beams and re-convolved with a Gaussian beam) with the 
Maximum Entropy Method (MEM), which is more suit- 
able for the extended signals than CLEAN (e.g., Narayan 
& Nityananda 1986; Sault 1990). Non-linear deconvolu- 
tion is applied using the MIRIAD task mosmem to the 
dirty images either separately or jointly with the single- 
dish image. Since the SZE brightness is negative at i^ < 220 
GHz while MEM is applicable only to the positive sig- 
nals, we change the sign of the entire maps when applying 
MEM. The pixel size is taken to be the same as the input 
image (1^^). 

Since a proper estimate of the underlying noise level of 
an image with an extended source is not straightforward, 
we adopt and compare the following two definitions us- 
ing the pixels within 120'' x 120'' (Models A and C) or 
240" X 240" (Model B); 1) ath is the rms value of a decon- 
volved image with the same observing parameters but us- 
ing the null input model including only the thermal noise, 
and 2) (Jdiff is the rms value of the difference map between 
the deconvolved output and the smoothed input; the in- 
put described in §2.1 is convolved by a Gaussian with the 
same FWHM as the synthesized beam of the output. The 
values of ath are given in Table 2. Thermal noise is nearly 
constant over the map except near the edge. 

Finally, we quantify the quality of reconstruction by the 
image fidelity defined by (Pety, Gueth & Guilloteau 2001) 



102 



f{0) 



|jsmooth(^)| 



max{|/out(^) - I!r"'''i^)l O.Tadiff} 



(5) 



where /out(^) is the output intensity of the deconvolved 
image at a sky position 6 and /fj^°°*^(^) is the input inten- 
sity smoothed with the same synthesized Gaussian beam 
as the output. The above fidelity roughly corresponds to 
an inverse of the fractional error of reconstruction. The 
second term in the denominator is introduced to avoid 
an overshoot of the fidelity due to a coincidental match 
between the input and the output. 

3. Results 

3.1. Model A: Gaussian 

We plot in Figure 3 the results of a representative case 
of a Gaussian model with ^modei = 40" FWHM after total 
integrations of 2 hr 12mx50, 7mxl2, and SDx4, respec- 
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Fig. 4. Mean fidelity of the deconvolved images as a function 
of ^model ^OT Modcl A. The mean is taken over the central 
70^^ X 70^^ and ^model x ^model for SDx4 and the others, re- 
spectively. Symbols indicate the results for 12m x 50 (crosses), 
7mxl2 (open triangles), SDx4 (open circles), and 12m x 50 + 
7mxl2 + SDx4 (filled circles). 



tively. With the 12mx50 data alone, reconstruction is 
poor and extended (^ arcsec) feature is largely missed 
even though the noise level of the 12mx50 image (panel 
b) is ~ 1/100 of the peak signal. This is also indicated 
by the fact that <Jdiff ^ crth- On the other hand, signif- 
icant improvement in the image fidelity is achieved once 
low spatial frequency data from 7mxl2 and SDx4 are in- 
cluded (panels d, e, f). The brightness of the SDx4 image 
appears to be low simply because it is diluted by a large 
beam (panel c). 

Figure 4 illustrates how the mean fidelity around the 
center varies with the spatial extent of the emission. It 
is evident that the interferometric data lose sensitivity 
for the emission much more extended than the size of 
their synthesized beams, i.e., ^beam ^ 5^' for 12m x 50 and 
^beam ^ 20^' for 7mxl2. The single-dish data play an es- 
sential role in improving the sensitivity at larger scales. 

3.2. Model B: Bullet cluster 

Figure 5 shows the images of a simulated bullet clus- 
ter after total integrations of 10 hr for 12mx50 and 40 hr 
for 7mxl2 and SDx4. We have fiducially fixed the ra- 
tio between the integration times at its nominal value of 
1:4 (Pety, Gueth & Guilloteau 2001; Morita & Holdaway 
2005). The impact of changing the ratio will be discussed 
in §3.5. 

The morphology of the input SZE image (Fig. 5a) is 
quite different from that in X-rays plotted in Figure 1 of 
Akahori & Yoshikawa (2012). In particular, the so-called 
"bullet" , the X-ray peak of the sub-cluster, is not apparent 
in the SZE image as it lies near the contact discontinuity 
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Fig. 6. Comparison of the input model (thick dashed hne) 
and the mock data (error bars) for Model B over a strip with 
width 10^^ shown in Figure 5(e). The bin size is 10^^ The 
vertical lines indicate the positions of the shock front (thin 
dashed) and the contact discontinuity (thin dotted), respec- 
tively, on the collision plane of the two merging clusters. 

across which pressure is nearly constant. On the other 
hand, the shock front, i.e., the pressure gap, ahead of 
the bullet is more prominent (see also Figs 6 and 8) and 
extends over several hundred kpc. The SZE peak lies near 
the X-ray peak of the main cluster where electron density 
and temperature are both high. 

There are 12 radio point sources over 3lC x 3lC among 
which 10 lie within the 240'^ x 240^' region mapped by 
12m X 50 (Fig. 5a). All the sources brighter than '-«^0.1 
mJy are detected above a conservative threshold of 7 ath 
in the 12m x 50 data as indicated in Table 3 and Figure 
5(b). After correcting for the side- lobe effects as described 
in §2.3, reconstructed fluxes agree with the input values 
within 20% for the resolved point sources and the source 
positions are identified more accurately. While the aper- 
ture bias of the source fluxes has been enhanced by point- 
ing errors, it is still comparable to the thermal noise and 
we do not attempt to correct it further in the present 
paper. Two bright sources in the guard band are also 
detected, although the noise level there is larger due to 
the primary beam attenuation beyond the map edge. The 
remaining 3 fainter sources are undetected but their con- 
tamination to the SZE is negligible. 

Once the detected sources are removed, the deconvolved 
image by 12mx50 + 7mxl2 + SDx4 reproduces the input 
model with the maximum fidelity reaching ~ 16 (panels 
e and f in Fig. 5). The fact that ath — crdiff holds also 
assures that the reconstruction of the extended signal is 
as good as expected. 

Figure 6 further illustrates one-dimensional SZE inten- 
sity profile across the shock front. The error bars indicate 



where a^td is the standard deviation of the pixel data in 
the bin and A/'beam = ^binMbeam is the number of syn- 
thesized beams contained in the bin area Abin- In equa- 
tion (6), we use A'beam, instead of the number of pixels, 
because the pixel values in a deconvolved image are cor- 
related over the beam area Abeam- Note also that abin 
does not include explicitly the thermal noise of the single- 
dish data (see Table 2) that is responsible for fluctuations 
on spatial scales larger than Abin- It is evident from this 
figure that the overall structure of this cluster is well re- 
constructed, although the shock front is only marginally 
resolved. We therefore explore the feasibility of a targeted 
observation toward the shock below. 

3.3. Model C: Shock front 

We plot in Figure 7 the images toward the shock front 
after total integrations of 10 hr for 12m x 50 and 40 hr for 
7mxl2 and SDx4. Again, reconstruction is significantly 
improved by adding the short baseline data of 7mx 12 and 
SDx4. 

Figure 8 shows that the high angular resolution SZE 
observation with ALMA is a powerful tool to resolve the 
shock front, characterized by the temperature and pres- 
sure jumps. On the other hand, the X-ray emission is weak 
in the low-density shock region and rises sharply behind 
the contact discontinuity. The Compton y-parameter or 
projected pressure is nearly constant behind the contact 
discontinuity (a slight disagreement between its peak and 
the contact discontinuity is due to the projection effect). 
The SZE and X-rays are thus complementary in probing 
the detailed shock structure and the former is particularly 
useful in detecting hot rarefied gas. 

The width of the shock heated region is ~ 60 kpc cor- 
responding to ~ 14^' dX z = 0.3 as displayed in Figure 8 
or ~ 8^' at z = 1. This physical size is consistent with 
the X-ray data of this cluster presented in Figure 32 of 
Markevitch & Vikhlinin (2007), whereas the errors in the 
inferred X-ray spectroscopic temperatures are still large. 
Given the fact that the intensity of the SZE is redshift 
independent and the angular diameter distance is nearly 
constant at z > 1, our results strongly point to the good 
capability of ALMA in resolving shocks at such high red- 
shifts. 

3.4' Systematic errors of the single- dish data 

As noted in §3.1, single-dish data play an essential role 
in reconstructing the extended signals. Due to a larger 
beam size and a smaller number of antennas, the single- 
dish data are likely to be affected more severely by sys- 
tematic errors, such as calibration errors and atmospheric 
fluctuations, than the interferometric data. To check the 
impacts of uniform and random systematic uncertainties, 
we re-run the mock observations of Model C by i) reducing 
the SDx4 gain by 20% or ii) enhancing the random noise 
of the SD X 4 data by an order of magnitude. For the latter. 



Yamada et al. 



[Vol. , 




20 -20 

RA offset [arcsec] 




20 -20 

RA offset [arcsec] 



20 



10 



t 

cd 



-10 



-20 



1 1 1 1 1 1 1 1 1 1 


1 1 i| 1 1 1 


1 - 


<n > x""^s 

<--'■ ■■■■/ 


\ 1 T 

\ 1 spec 


— 


- Ix 


\ 1 

N 1 

N 1 

N 1 

N 

\ 1 

] 

"-4 


---: 


/ 

/i 

/ 1 

/ 1 

/' 1 

/ 1 


- 


-y 


/ 1 


; 


1 1 1 1 1 1 1 1 1 1 


1 III 1 1 1 


1 




40 20 -20 

RA offset [arcsec] 



-40 



40 20 -20 -40 

RA offset [arcsec] 



Fig. 8. Top: Same as Fig. 6 except that the results of Model 
C shown in Fig. 7(e) are plotted with a smaller bin size of 
5^^ Bottom: Physical quantities with arbitrary units of the 
input model for the same region as plotted in the top panel. 
Thick lines indicate the projected temperature proposed by 
Mazzotta et al. (2004) that mimics the X-ray spectroscopic 
temperature (solid), the mass weighted electron density (short 
dashed), the thermal bremsstrahlung X-ray brightness (long 
dashed), and the Compton y-parameter with a negative sign 
added for comparison with the SZE decrement (dot-dashed). 



we also examine a linear method, often called feathering, 
to combine the single-dish data instead of non-linear joint 
deconvolution; a deconvolved image is produced only for 
12mx50+7mxl2 by MEM and then combined with the 
SDx4 image using the MIRIAD task immerge. 

Figure 9 shows that introducing a reduction in the 
single-dish gain by ~ 20% results in a systematic decrease 



Fig. 9. Same as the top panel of Fig. 8 except for introducing 
a systematic reduction in the single-dish gain by 20% (top) 
or additional random noise on the single-dish data that is an 
order of magnitude larger than the thermal noise (bottom). 
Triangles in the bottom panel show the results of the lin- 
ear method of combining the single-dish data (error bars are 
omitted for clarity). Error bars do not explicitly include the 
single-dish noise as described in the text. 

in the reconstructed brightness by ~ 20% as well. On the 
other hand, the enhanced random noise primarily leads to 
an offset in the zero brightness level, whereas the shape 
of the emission profile is nearly unchanged. The latter 
bias has been known as an intrinsic problem with MEM 
for low S/N data (e.g., Cornwell, Braun & Briggs 1999). 
The linear method described above gives a better estima- 
tion of the zero level but poorer results on the emission 
profile (bottom panel). This is because the linear method 
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Fig. 10. Mean fidelity of the mock bullet cluster (Model B 
without point sources) as a function of the integration time of 
AC A (7mxl2 and SDx4). The integration time of 12m x 50 is 
fixed at 10 hr. The mean is taken within the central 60^^ x 60^^ 
in the 12m x 50 + 7mxl2 + SDx4 image without smoothing 
(circles) or after being smoothed to the effective beam FWHM 
of 10^^ (crosses) and 30^^ (triangles). 

directly combines the single-dish data with the interfer- 
ometer data on the u — v plane and hence retains noisy 
pixel values in the single-dish map. On the other hand, 
MEM searches for a smooth solution that accounts for the 
global structure of the single-dish data and is less sensitive 
to individual pixel values. In either case, the above results 
suggest that the quality of the SDx4 data may limit the 
accuracy of the SZE observations. 

3.5. Relative integration time between the 12m array and 
ACA 

In Models B and C, the ratio between the integration 
times of 12mx50 and ACA (7mxl2 and SDx4) has been 
taken at its nominal value of tACA/^i2m = 4 (Pety, Gueth 
& Guilloteau 2001; Morita & Holdaway 2005). We now 
examine the impacts of varying this ratio by re-running 
wide-field Model B simulations. The total integration 
time of 12mx50 is fixed at 10 hr and the point sources 
are excluded from the input for simplicity. 

Figure 10 indicates that tACA/^i2m ^ 1 is required to 
achieve a ~ 10% accuracy in the image reconstruction 
at the resolution of 12mx50, Oheam ~ 5'^ For smaller 
^ACA/^i2m, the thermal noise of ACA in /iJy/arcsec^ dom- 
inates over that of 12mx50, as can be readily expected 
from Table 2. For larger tACA/^i2m, on the other hand, 
the image fidelity increases more rapidly at larger scales. 
The optimal value of tACA/^i2m thus depends on the spa- 
tial scale of interest, and the nominal value of 4 appears 
to be a reasonable choice for the case considered in the 
present paper. 



4. Discussion 

High resolution SZE mapping with ALMA has pro- 
found implications on cosmology and structure formation. 
First of all, the shock velocity inferred via the Rankine- 
Hugoniot jump condition from Chandra for IE 0657-558 
amounts to --4700 km s"^ (Markevitch & Vikhlinin 2007), 
which appears to be exceptionally high within the frame- 
work of a concordance ACDM model (e.g., Hayashi & 
White 2006; Lee & Komatsu 2010). Our results suggest 
that ALMA can significantly improve both quality and 
quantity of such studies by detecting shocks up to z > 1 
and measuring directly the pressure gap across the shock. 

The existence of very hot electrons with /cTe ^10 keV 
in clusters has also been inferred by the past SZE obser- 
vation (Kitayama et al. 2004) and the hard X-ray spec- 
troscopy (Ota et al. 2008; Nakazawa et al. 2009; Sugawara, 
Takizawa & Nakazawa 2009). Nevertheless their spatial 
distribution as well as the link to the non-thermal compo- 
nent are still unclear. ALMA will be highly complemen- 
tary in probing the nature of these relativistic electrons to 
hard X-ray missions such as NuSTAR^ and ASTRO-H^, 
whose spatial resolutions are 45^' and 1.7^ respectively. 

Our results can further be used to estimate the feasi- 
bility of SZE observations in other bands. Frequency de- 
pendence of the SZE maps will provide a useful probe of 
line-of-sight gas motions via the kinematic SZE, electron 
temperature via relativistic corrections, and the existence 
of non-thermal electrons (e.g., Colafrancesco, Marchegiani 
& Buonanno 2011; Prokhorov et al. 2011). The SZE 
increment at frequencies above 600 GHz has also been 
detected by Herschel toward 1E0657-558 (Zemcov et al. 
2010). 

For given arrays and bandwidths, the thermal noise in 
Jy/sr is proportional to T^yst^^^ v^ ., where Tgys is the sys- 
tem temperature, tint is the integration time and v is the 
observing frequency. If one is to observe a source of inten- 
sity /src in Jy/sr with A^mos mosaics and perform smooth- 
ing over A^sm beams, the signal-to- noise ratio is given by 



S_ 

N 



(X 



^srcl^^int 



1/2 ^1/2 



f- sys^ 



N, 



1/2- 



(7) 



Given that both A^mos and A^sm are nearly proportional to 
u'^ for fixed observing area and effective beam size after 
smoothing, the integration time to reach a given signal- 
to-noise ratio is 



tu 






(8) 



Figure 11 shows the quantity given in equation (8) in 
the cases of the thermal SZE with /cTe = 15 keV and the 
kinematic SZE with the line of sight velocity of F = 3000 
km s~^, including relativistic corrections (Itoh & Nozawa 
2004; Nozawa, Itoh & Kohyama 2005). The vertical axis 
shows the ratio with respective to the integration time for 
the thermal SZE at 90GHz. We use system temperatures 



^ http://www.nustar.caltech.edu/ 
^ http://astro-h.isas.jaxa.jp/ 
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Fig. 11. The ratio of integration times, as a function of ob- 
serving frequency, required to reach a given signal-to-noise 
ratio in a fixed effective spatial resolution over a fixed ob- 
serving area. Symbols indicate the cases for the thermal SZE 
with kTe = 15 keV (triangles) and the kinematic SZE with 
V = 3000 km s"-*^ (crosses). The vertical axis is normalized 
with respective to the integration time for the thermal SZE 
at 90GHz. 

at zenith (RA,DEC) = (00 : 00 : 00, -23 : 00 : 00) given 
by the ALMA sensitivity calculator. Precipitable water 
vapour is taken to be 2.748 mm at i^ < 163 GHz, 1.796 
mm at 163 < z/ < 275 GHz, and 1.262 mm at 27b<iy< 500 
GHz. The resulting system temperatures tend to increase 

^40GHztoTsys-180Kat 



gradually from Tgys '^ 40 K at 



u ^ 350 GHz, with some overshoots wherever atmospheric 
opacity gets high, e.g., at z/ ^ 180 GHz. The rise of the 
integration time for the thermal SZE at z^ ~ 220 GHz is 
due to the null of its intrinsic spectrum. It is obvious that 
the SZE observations, particularly that of the kinematic 
SZE, become highly time-consuming at z^ > lOOGHz due 
both to an increasing Tgys and to the z^^ factor in equation 
(8). On the other hand, much faster mapping is possible 
at lower frequencies, e.g., z^ = 40 GHz. 

Finally, we estimate the impact of unresolved (< 0.1 
mJy) point-like sources as follows. The 150GHz source 
counts from recent SPT and ACT observations are ap- 
proximately A^(> Su) ~ 2 X 10~^(5'i5o/10mJy)~^ arcmin"^ 
for synchrotron-dominated sources and N{> Sj^) ~ 
10~^(5'i5o/10 mJy)~^-^ arcmin"^ for dust-dominated 
sources at 5'i5o ~ 10 mJy (Vieira et al. 2010; Marriage 
et al. 2011). If we simply extrapolate these relations to 
fainter fluxes at 90 GHz, the nominal confusion criterion 
of 1/30 beam~^ , or equivalently 4 arcmin"^ for 6>beam = 5'^ 
is reached at Sqq ~ 0.6 /iJy and 10 /iJy, assuming the av- 
erage spectral shapes of Sj^ ex u~^-^ and u^ for the two 
populations (e.g., Vieira et al. 2010), respectively. This 
suggests that the source confusion is not likely to affect 
severely the SZE observations considered in this paper. 



unless there is a significant excess of faint sources with 
respect to the above extrapolation. 

5. Conclusions 

We have performed imaging simulations of the SZE of 
galaxy clusters with ALMA including ACA. In its most 
compact configuration at 90GHz, ALMA will resolve the 
intracluster medium with an effective angular resolution 
of ~ 5'^ corresponding to ~ 20 kpc at 2: = 0.3 and ~ 40 kpc 
at z = 1. Such observations will be particularly useful in 
detecting shock fronts and/or relativistic (either thermal 
with kTe ^ 10 keV or non-thermal) electrons produced 
during violent cluster mergers at high redshifts, that are 
hard to resolve by current and near-future X-ray detec- 
tors. 

Our results imply that ACA plays an essential role in re- 
constructing the global structures of the SZE and its capa- 
bility may be limited by the accuracy of single-dish data. 
Expected sensitivity of the 12m array based on the ther- 
mal noise is not valid for the SZE mapping unless accom- 
panied by an ACA observation of at least equal duration. 
An optimal ratio of integration times between the two ob- 
servations depends on the spatial scale of interest and the 
nominal value of 4 appears to be a reasonable choice in 
the present case. The SZE observations at z^ > lOOGHz 
will become excessively time-consuming owing to the nar- 
rower beam size and the larger system temperature. On 
the other hand, significant improvement of the observing 
efficiency is expected once Band 1 is implemented in the 
future. 

In addition to the SZE, ALMA wih detect radio/IR 
sources in the field-of-view. The spectral energy distribu- 
tions of such sources in the ALMA bands are still highly 
uncertain and will be crucial both in understanding the 
nature of the sources and in quantifying their contamina- 
tion to the SZE survey data (e.g., Lin et al. 2009). 
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Note Added in Proof (June 22, 2012) 

After this paper was accepted for publication, updates 
to the ALMA sensitivity calculator released on May 31, 
2012 resulted in an increase in the system temperature at 
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90GHz by 10% and an inclusion of quantization efficiency 
of 0.96, compared to those listed in Table 1. These mod- 
ifications would simply raise the thermal noise amplitude 
by 14% and the results of the present paper remain es- 
sentially unchanged once the quoted integration times are 
scaled upward by 30%. 

Appendix 1. Effective Integration Time in 
Mosaicing Observations 

In mosaicing observations, effective integration time to- 
ward a given sky position 9 is expressed as 

^int 



ieff (^) = ei,{e) 



N„ 



(Al) 



where tint is the total integration time over the entire map- 
ping area, A^mos is the number of mosaics, and €^{6) is 
the "mosaicing overlap factor" introduced in Holdaway & 
Rupen (1995) and Morita & Holdaway (2005). The last 
quantity accounts for the overlap of mosaics that cover 
the same position and is given by 






eiie)=yg\0-0p), 



(A2) 



where g is the primary beam gain and Op is the pointing 
center of the p-th mosaic. In what follows, we derive the 
values of (e^) and (teff) analytically. 

First, we consider primary beam gain functions of the 
form 



g{0) 



exp 



V 2^y 



2Ji(7rD6>/A) 
-kDO/X 



Gaussian, 



Airy disk. 



(A3) 



where Ji is the Bessel function of the the first kind of 
order unity, D is the diameter of the telescope, and A 
is the observing wavelength. Integrating ^^ over the sky 
gives the effective area of integration per pointing as 



2'A^/a 
7rD2 



: 0.567 ^2^ 



Gaussian, 



'. 0.554 $1^ Airy disk. 



where 



I A 



f 

Jo 



■MxY 



dx ci 0.0575, 



(A4) 



(A5) 



and ^pb is the primary beam FWHM corresponding to 
v81n2(7G for Gaussian and 1.028 A/I) for the Airy disk. 
For each p in equation (A2), g^ ^ I within Again centered 
at Op, whereas ^^ ^ otherwise. 

Second, we suppose that mosaicing centers are placed 
on a regular grid, either square or triangular, with con- 
stant spacing d. The physical area around each grid point 
enclosed by perpendicular bisectors with the neighboring 
points is 

r d'^ square grid. 



^d' 



where the shape of the enclosed area is a square and a 
hexagon, respectively. 

Finally, we average equation (A2) around an arbitrary 
grid point Oq that lies sufficiently far (^ \/ Again /tt) from 
the map edge, to obtain 



{< 



Ep=rI,9Ho-o,)d'e 



Igrid 



(A7) 



where J. denotes an integral around the i-th grid point 
over the area Ag^id- The last result in equation (A7) does 
not depend on q and remains unchanged even if the av- 
erage is taken over multiple grids. Substituting equations 
(A4) and (A6) into (A7) yields 

where the representative values in the case of Nyquist 
spacing are 



(^m)Ny 



2.27 Gaussian, square grid, 

2.62 Gaussian, triangular grid, 

2.22 Airy disk, square grid, 

2.56 Airy disk, triangular grid. 



Equations (Al) and (A7) give 

(^eff) (Cm) ^gain 



^int -^^mos ^map 



(A9) 



(AlO) 



where A^mosAgrid roughly equals the mapping area Amap- 
Equation (AlO) indicates that the mean integration time 
at each sky position within fixed A^ap does not (apart 
from the time spent at the map edge) depend on d or 
the grid orientation; (e^) and Amos are both proportional 
to d~'^. On the other hand, its dispersion does as shown 
below. 

We further compute numerically the mean and the 
standard deviation of e^ by creating square and trian- 
gular grids over the sky area with 56>pb x 56>pb and ex- 
tracting 10000 random positions of from the central 
3^pb X 3^pb. Figure 12 illustrates that Ae^/(e^), or equiv- 
alently Ateff/(teff), drops sharply to <C 1% at d < 0.5^pb 
in all the cases considered here. We have checked that 
(^m) computed numerically agrees with equations (A8) 
and (A9) to better than 1% and 4% at c^ < 0.76>pb and 
d < 6>pb, respectively. The agreement improves further 
once the larger computational area is adopted. 

The above results imply that, as far as the Nyquist con- 
dition is satisfied, square and triangular grids lead to prac- 
tically similar and uniform effective integration time apart 
from the map edge. The difference due to assumed shapes 
of the primary beam, Gaussian or an Airy disk, is not sig- 
nificant either. 
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axes are normalized by the mean (e^) and the primary beam 
FWHM ^pb, respectively. 
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Fig. 2. The u — v coverages of our mock observations with 12mx50 (left) and 7mxl2 (right) of Models A (top), B (middle), and 
C (bottom). 
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Fig. 3. Mock images of Model A (Gaussian) with ^model = 40^^ Except for panel (f), the color scale is in /x Jy/arcsec^ . A hatched 
ellipse at the bottom-right corner indicates the synthesized beam, (a) Input model convolved with the same synthesized beam as 
panel (e). (b) Deconvolved image for 12mx50 with (crth,crdiff) = (1.1,17) [//Jy /arcsec^] and ^beam = 3.9". (c) SDx4 image with 
<^th = 0-64 /iJy/arcsec-^ and ^beam = 69^^ (d) Deconvolved image for 7mxl2 + SDx4 with (crthjCr^iff) = (0.48,0.79) [/iJy /arcsec-^] 
and ^beam = ^^" • (e) Deconvolved image for 12mx50 + 7mxl2 + SDx4 with (ath,f^diff) = (1.0,1.2) [/iJy/arcsec^] and ^beam = 4.0^^ 
(f) Fidelity of the image shown in panel (e). 
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Fig. 5. Similar to Fig. 3 but for Model B (bullet cluster), (a) Input model convolved with the same synthesized beam as panel 
(e). Symbols indicate the positions of the X-ray peaks of the main cluster (plus) and the sub-cluster (cross), and the shock front 
(square), (b) Dirty image of 12m x 50 with crth = 0.62 /iJy/arcsec^ and ^beam = 4.5^^ Circles indicate the point sources that are 
detected and removed in the subsequent panels, whereas squares denote undetected ones, (c) SDx4 image with crth = 0.23 /^Jy/arcsec-^ 
and ^beam = 69^^ (d) Deconvolved image for 7mxl2 + SDx4 with (crth,crdiff) = (0.15,0.18) [/xJy /arcsec^] and ^beam = 19^^- (e) 
Deconvolved image for 12mx50 + 7mxl2 + SDx4 with (crth,crdiff) = (0.52,0.53) [/iJy/arcsec^] and ^beam = 4.9^^ A box indicates 
the region over which the profile in Fig. 6 is computed, (f) Fidelity of the image shown in panel (e). 
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Fig. 7. Similar to Fig. 3 but for Model C (shock front), (a) Input model convolved with the same synthesized beam as panel 
(e). (b) Deconvolved image for 12mx50 with (crth, CTdiff), = (0.35,0.95) [/iJy /arcsec^] and ^beam = 4.5^^ (c) SDx4 image with 
(7th =0.14 /iJy/arcsec^ and ^beam = 69'^ (d) Deconvolved image for 7mxl2 + SDx4 with ((Tth,(^diff),= (0.099,0.18) [/iJy /arcsec^] 
and ^beam = 18^^- (e) Deconvolved image for 12m x 50 + 7mxl2 + SDx4 with (crth,crdiff) = (0.30,0.33) [/iJy/arcsec^] and ^beam=4.8^^ 
A box indicates the region over which the profile in Fig. 8 is computed, (f) Fidelity of the image shown in panel (e). 



